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THEORETICAL STUDIES ON THE MECHANICAL BEHAVIOR OF 
GRANULAR MATERIALS UNDER VERY LOW INTERGRANULAR STRESSES 


by 


Kenneth W. French, Jr. 
Professor of Mechanical Engineering 
John Brown University 
Siloam Springs, Arkansas 


ABSTRACT 


The report describes the salient aspects of the theoretical modeling of a 
conventional triaxial test (CTC) of a cohesionless granular medium with stress 
and strain rate loading. Included are a controllable gravitational body force 
and provision for low confining pressure and/or very low intergranular stress. 
The modeling Includes rational, analytic, and numerical phases, all in various 
stages of development. 

This work is an extension of two previous summer’s contributions to an ongoing 
project at NASA/MSFC and with contractors at the University of Colorado. Hence, 
this report builds strongly on the base of those from 1984 and 1985, The 
experimental culmination of the project will be performed in a mid-deck locker 
on a future Space Shuttle Flight. The numerical evolutions of theoretical 
models will be used in final design stages and in the analysis of the 
experimental data. 

In this the experimental design stage, it is of special interest to include in 
the candidate considerations every anomaly found in preliminary terrestrial 
experimentation. Most of the anomalies will be eliminated by design or 
enhanced for measurement as the project progresses. The main aspect of design 
being not the physical apparatus but the type and trajectories of loading 
elected. The major considerations that have been treated are; appearance and 
growth of local surface aberrations, stress-power coefficients, strain types, 
optical strain, radial bead migration, and measures of rotation for the proper 
stress flux. 
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INTRODUCTION 




The 1980s seem to be the deuiiue for the coming cf cgs for advanced enelysls In 
rheology and soil mechanics. They are quickly catching up to their simpler, 
faster developing sisters; fluid and structural mechanics. In the granular 
mechanics community there seems to be almost a rush to install the full results 
of rational modeling and large strain plasticity theory. Rates of stress and 
strain are being incorporated into constitutions, objectivity is demanded, 
stress power is accounted, and damage, memory, and fracture are actively 
courted. However, there is something grander about this whole development. 

The behavior of granular material has long escaped analytic capture for a very 
good reason; it exhibits a very complicated constitution. Granular materials 
have evaded analysis except in the very restricted local sense of using test 
coefficients in an interpolated analytic framework. Often even then a 
confining history must be specified. Further the dissimilarities among soil 
types is wondrously wide, to say nothing of the spectrum of man-made substances 
with discrete particles from elastic /brittle glass to polymeric types with 
memory t 


The new models are immediately followed by computer analysis which by 
prediction enhances final design stages and which by correlation enhances the 
results of prototype and final experiments. Nevertheless as the lessons in 
turbulence in fluid mechanics and buckling in structural mechanics have shown, 
the progress with numerical analysis will still prove to be a frontier for many 
years • 

Whether the bevy of classical soil mechanics tests are thought of as studies of 
the behavior of infinite half-spaces of granular material or as a search of 
constitution by partial derivative or as preparation for the use of membrances 
for control of granular bulk; the influence of gravity body force is crucial to 
a full understanding. Our human intuition is built on experiences with 
granular media under the domination of gravity . Outside that realm there are 
precious few (albeit extremely interesting and crucial) cases and those highly 
transient. Separation by flow of multicomponent saturated river beds, 
liquefaction in quicksand or under earthquake or other strong cyclic loading, 
cyclic mobility deformations; but at best these are transient snapshots of low 
intergranular stress behavior. 

There is also the very interesting prospect of the granule being the ideal form 
for storage of bulk material in zero-g. The underlying reason being its ability 
to exhibit high or low frictional forces. Unfortunately, this characteristic 
is strongly inhibited by strong gravitational forces which confine the granules; 
hence, soil liquefaction is a rarity on earth which does not sadden civil 
engineers . 
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Vibratory transport with fluxes of discrete masses has been 
employed commercially for many years. This in combination with an 
oscillating confining pressure has potential; the former to move the 
latter to act as a check valve. In space a double membrane container 
with pressure control of the intermembrane region and vacuum or 
pressure control to the inner region would allow behavoir rangeing 
from near fluid to tremendously frictional and hence either pumping 
or the retention of the integrity of an almost arbitrary shape. 

Therefore either bulk storage or peristaltic transport are realizable 
within the same containment. That is a granular material does not 
have the high bonding of a frozen solid nor is it dominated by surface 
tension or vapor pressure typical to liquids - 

The concept of a granular media as a 3-D fabric with noncatastrophic 
damage being done continuously by the loading represents perhaps 
the best thermodynamic view and keys on vision of stress-power. 

Ultimately this damage becomes particle crushing, shear banding, fabric 
fracture and if loading is strongly cyclic or containment small, cyclic 
mobility or liquif action. 

It is into this fascinating realm that the "Costes Team” 
commenced finding a niche in 1978. The particular geometry chosen is 
the classical triaxial test with provision for saturation draining 
and a variety of loading trajectories; see Figures 1 and 2. The 
following section will detail the exact position of this work in the 
general development of the ongoing project. The experimental culmination 
being Space Shuttle borne tests of a wide variety and the constitutive 
theory an unending frontier of discovery based on these better low 
intergranular stress results. 



TRIAXIAL TEST 

FIGURE 1 FIGURE 2 
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OBJECTIVES 


The primary goal of this work is to provide theoretical 
modeling support for the Mechanics of Granular Media ( MGM ) , 
activity with special direction toward the models for the 
conventional triaxial configuration with specimens of cohesionless 
spheroidal material. In this respect it is a true continuation 
of the past years' work. 

The objectives of this work are to provide state of the art 
rational modeling, analysis within the precision required, and 
numerical support for the full range of constitutions and for 
anomalies outside the model. Typically the model is local, 
a differential continuum; the analysis is global, integral; and 
the numerical approach attempts both local and global approximations. 
It is necessary to explicitly include acceleration of gravity and 
mass density. The former for body force control and the latter 
for that as well as density sensitivity in the constitution and 
for tracking dilation. 

One of the primary objectives is to confront the real boundary 
conditions of platen encroachment, membrane adjustment and 
penetration, and boundary shear stress power and to resign them 
to the constructions used in analysis. This is crucial; to relate 
all the statements of deformation and its time and space derivatives 
through all three theoretical phases. The careful, definitive 
statement of strain will permit a correct extension to the use of 
optically measured displacements on the cylindrical boundary. 
Preliminary calculations of the inherent experimental strain 
uncertainty are desireable to anticipate the necessary prediction 
level for the analytic and numerical phases. 

The numerical objective is to provide a local version of the 
MICROFEM code that is fully interactive, specialized to a reasonable 
level of entry data, implemented with screen prompting and 
segregated from the FEM skeketon. A brief manual will be produced 
to guide operation and to document the flow of the code through 
subroutines so that modification can be made and behavoir traced to 
its origin. This code should accomodate a variety of constitutive 
models, but be limited to a single solver and to the geometry, 
and loading trajectories envisioned for the culmination experiments. 
The FEM must reenter the constitutive model at every iteration and 
be fully sensitive to the boundaries of behavoir internal to the 
stress space. 
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JAUMANN QUANDARY 


A conference paper in 1981 showed that under certain 
conditions the combination of a Prager-Ziegler kinematic hardening 
rule with a Jaximann stress flux can cause an oscillatory stress- 
strain relationship. Figures 3 & 4 show some extreme and more mild 
forms of this oscillation; they are from Lee (82) and Dafalias(82) 
respectively. It is interesting to note that Lee et al were the 
editors of the conference at which the original paper was presented. 

The conclusion of Lee (82) and Dafalias(82) and subsequently 
of Kiousis(86) was that the Jaumann stress flux was inappropriate in 
that particular circumstance and others. Kiousis states that the 
Jaumann stress "rate” ( stress flux) is not applicable to materials 
exhibiting kinematic hardening. Whether this conclusion is verified 
or if the constitutive model assumed is at fault remains to be seen. 

Since there is a wide variety of choices for the stress flux, 
it is omnious to presume that one must be selected for each combination 
of constitution and loading trajectory. Eringen(80) offers seven 
permutations. It is noteible that the property of objectivity is 
distributive to transformation so new versions can be formed by adding 


objective 

tensors or groups. 
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The basic set of stress flxaxes are built from grad V , and stem 
from the relation between material and moving coordinate frames. 

The need arises when incorporating stress rate ( material derivative) 
into the constitutive relation because stress rate is not inherently 
objective. This usually amounts to a sort of Boussinesq approximation 
in that grad V is ignored in the momentum equation but included in 
the constitutive statement. Thus the momentiam equation becomes: 
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In order to lay a groundwork for the new stress fluxes, the 
following is a condensation of Truesdell (77) II chapters 5-11. 


v\/ (G = /D ^ W 


velocity gradient 
deformation gradient 


polar decomposition . . rotation tensor , RH stretch tensot 

LH stretch tensor 


(E* 


< V is frame indifferent ^ 

V = JR1LJP(^ 

^ ILT = JetV 

LT* = * "V** FF^ 


F ht4jR 


are not 


jacobian 


carries principal axes of strain in X 
into principal axes of strain in x 


I ; TLT=V ^ 


pure £sti. etch 


W • * iR (uij‘- 

Spin components from polar decomposition 


ID ^ iJRCURJ t 

deformation rate components from polar decomposition 
orthogonal tensor between frames 

£>* " (8 £) (Q (hence objective) 

b«t= w* = d)W A. 


A 


is the spin of frames joined by 
relative to one another 

so that the asterik deformation rate is frame indifferent vrtiile 
the asterik spin is the original plus the relative. 
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Isotropic 

hardening 
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FIGURE 3 LEE STRESS FLUX FIGURE 4 DAFALIAS STRESS FLUX 




Considering a fairly general constitution: ( eg Scott (84)) 

«• firfl-ine = Cjp *C Jjy *C -3^1 

i? ijUA £>XUlli * j. 

Entering the arena with just the first term on the constitutive 
side, Lee (82) and Dafalias(82) independently sought an acceptcd)le 
candidate for SI or beyond to a flux side that does not contain 
Jaumann terms explicitly. Their suggestions ( see Figures 3 & 4 ) 
resolve the problem in the sense that the oscillations are suppressed 
for the constitution and loading they considered: 


n 

JRJFi' 

. . . Dafalias(82) 

xn -» 

■w* 

. . . Lee (82) 


While this is satisfying in the immediate sense; it seems to 
lack the uniqueness vdiich is critical to any long range modeling 
progress. In fact it is in the worst spirit in the rational 
mechanics tradition continue such a strategy: 

1. preconceive a notion of acceptable results 

2 . select a reduced constitution 

3. "discover" a stress flux to satisfy #1. 

The use of such a strategy is very successful in early stages of 
discovery, in fact it is mode at which humans excel. However it 
is strictly for interpolation and often only for the physical and 
time scales for which the human has intuition. 

In spite of the apparent failure of rational mechanics to 
supply any direction, the most suspect portion of the model is the 
constitutive side. Following Einstein. .. "reduce to the simplest 
model but no further" .. .it seems that the constitution is in a 
state of over-reduction. In the interim this type of modeling; 

* d* ID 

has fluorished with the more progressive researchers in granular 
mechanics. Apparently none finding or at least reporting any 
ranges of oscillation or other nonphysical behavoir . On the other 
hand it is not difficult to find experimental results with what 
might be oscillation in the large strain region. Prevost(81). 

But of course the uncertainty in that realm masks such an observation 
from acceptance. 

While objective stress fluxes are legion, the strategy of 
selecting a constitution then shopping for a form of the stress flux 
that exhibits expected behavoir is a hedge against using the rational 
approach at all. Perhaps a thermodynamic ruling will have the same 
effect, but it must be explicitly discovered. 
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STRESS POWER 


In a general thermodynamic description of a material stress 
power ( nominally ) and mass density are related by the 

temperature variable between the "energy” and the "state" 
equations. The granular formulations while sometimes alluding to 
"adiabatic*,* invariably decouple the equations by removing any 
temperature sensitivity from the equations for energy and state. 

They are thus left related only by: 

e^=" jr * tr ® 

Density is occasionaly coupled to pressure ( 1/3 tr €T ) but the 
comparison is invariably on an integral, time independent basis 
( SR vs e ) - 

V 

Stress power and density are allied in that they appear in the 
coefficients of the full constitutive model.. [ ]... .Scott (84) 

following the general 10 coefficient model of Truesdell ( 55 ) ; 
basically this is a polynomial formulation in ^ and ID linearized in 
<T and ID alone but including the components of OID. There is perhaps 
a problem with this formulation for certain finite strain situations; 
consider the development below and notice that the Jaumann Quandary 
also becomes embroiled here for integrals over strain. 


TP - [ energy /voliime/ time] 

(ct: 'tJu -V- in ) nanin 

Jz I stres 

p • <r-H))r - ff-lD). V £ 

<cr - 


nominal 
stress power 1 


stress flux see J.Q. 


L J,3T *L * C 


constitutive 


Source Power: 


r*v 


Interstitial + Gravity+ Confinmentt Sink 

Jpj^ + Mji * * S 


kinetic energies being explicitly neglected. 


power sink terms both potential and dissipative 
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Power Sink Candidates: 


[ 1 ] 

[ 2 ] 


Platen Friction. . . 


deformation 


heat 


^def 


ritsiiujxcixie . 


• ^ luijd J. 


, / heat 

\ 


def 


hoop ^ 

heat 
^ def 

.penetration ^ 

^ heat 


[3] 


Fluid 



Pumping 

Compression 


Beyond the normal component for the encroaching platen which 
is the ultimate power source, the boundary working rate includes 
the following^ pumping and compression for the confining fluid, 
membrane stretching and penetration, and platen shear working rate. 
There is also the compression and pumping losses in the pore fluid, 
movement of the center of gravity ( if g ^ 0 ) and last but certainly 
not least the interstitial working rate of the granular medium. 

Another view of stress power is that commonly taken by the 
fracture mechanists, which is to assess damage as a quantitative 
measure of the "sinked" power; (ie the non elastic effects ..both 
the membrane and undrained fluids could also exhibit reversible 
compression ) and which appears in various forms in the coefficients 
of the constitutive relation. The fabric of the material is damaged 
in a way not delt with specifically , but its constitution changes 
in proportion to the s inked power. Kachanov (86) The use is sometimes 
made of a quantity termed, creep dissipation density; ^ ^ 

In the most general approach to this sort of analysis, Dunn & 
Serrin (85) introduce the interstitial working rate as: 

Both the creep dissipation and the interstitial working open the 
Jaumann Quandary once again. 
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STRAIN DEFINITION AND RESOLUTION 


This project has as a bounding condition a steadily 
encroaching platen with a centering force provided by a porus 
stone centrally located and of relatively small radius. The 
platen is nominally frictionless, hence at the moving platen 

boundary; . . , ^ . 

e = C, e = e^^ = full slip . 
zz rr @@ ^ 

That is the rates are prescribed independently of the rate 
behavoir brought into the motion by the constitution. Ultimately 
it is envisioned that the surface response will be sought as a 
sequence of optical images; optical segments will be identified and 
tracked. Particularly the images formed from the outer surface 
of the cylindrical surface. Subsequently an association will be 
made with the strain field and hence with the analytically 
generated strain predictions from the integrations of a continuum 
model . 


The following discussion is a consideration of the implications 
to the theoretical modeling of the anticipated order of magnitude 
of the movement and necessary transformations to establish a relation 
with optical displacements. This is an important precursor to the 
detailing of the constitutive model and its numerical counterpart. 

The intent is to identify and highlight the variety of ” strains** 
without resolving them all. This search is the prolegomenon to 
relating optical displacement, optical strain, computed strain, and 
true strain. That is the initiation of the correct processing 
algorithm and calibration requirements. 

Rates enter the specification from both the strobe nature of 
the optical technique and from the material constitution and loading 
trajectory (viz: constant strain rate ) ; the proper phasing of these 
three should be included in the design of the final experiment. 
Rotation of material lines, spin of embedded granules as associated 
with *' bulging and lumping **(see Thin Shell section ) is crucial 
as they will add to the distortion and must be compensated or 
the optical scheme must disallow these measurements automatically. 

As is discussed in the next section, the constitutive theory will 
not model the local surface phenomena, lumps, and probably will not 
predict the global surface phenomena, bulges. 

Some learning about evolution from the prototype experiments, 
may provide the prophylactic scheme for rejection of data during a 
surface abberation. For example; apparently local failures preceed 
global failure, and at the least their growth is enhanced by global 
failure. Hence the growth of local failure could predict location 
and age of global failure. 
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In Terzaghi & Peck (60) (ppl09) after giving a lengthy 
exposition the author adds "... from slow tests very much smaller 
values are obtained." In any venture of this sort the rate terms 
from the full constitution should be considered, if in no more 
detail than as order of magnitude studies of dimensionless weightings. 
Even the integral strain rate varies from its nominal value; 

e * C Vj) 

and the time rate for even increments of strain in that same sense; 

. / I “feVt ^ * 

At = (I- IT/ c-i*) 

It is also necessary to distinguish between a change in the time 
scales of response and true creep, because the local strain histories 
are quite different after appreciable strain. 
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Collecting the notation from Eringen(80); 


Lagrangian strain — 

1— 1 
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Eulerian strain- - 
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+ i (S’ 

mk 
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deformation gradients 

— 


*v 

= IE + IR 

= 

JF 




<C7u 

= IS + r 




velocity gradient - - 


Vv = 

= ID + W = 

= deformation 

rate 


with O indicating infinitesimal queintities of strain and rotation 


Hence vanishing infinitesimal strain are not sufficient to specify 
rigid deformation; infinitesimal strain is not a strain measure except 
under "infinitesimal" conditions, for vAiich Eulerian and Lagrangian 
strains are identical. 

• - ae 

fi " 5t & ® = d? + 


E = o infinitesimal deformations 

e = D if W<0(D) 

I e = P - eTV + We | 


it is also intereting that J - 1 + tr e 
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THIN SHELL CONSIDERATIONS 


There is an interesting possibility that a triaxial sample with 
a highly penetrated membrane exhibits a thin shell response for 
certain ranges of loading. The proximity of a tension element 
(membrane) and a compression element (granules) linked with inter- 
granular and membrane penetration friction forces sets the scene for 
thin shell behavoir. ( see Figure 5 ) Specifically the interest 
is only in the onset of buckling and post-buckling regimes. 

The interest in making these considerations arises for two 
reasons. First because of the appearance in testing of 3-D surface 
patterns which persist and grow in a nonhomogeneous way. Secondly 
because it is crucial to the rational model that in the region where 
theoretical model meets measurement there is a special physics. High 
penetration, large friction coefficient and their exponential coupling 
via membrane tension will allow the membrane to dominate the stress 
power, locally in a layer 0(d) thick. The suggestion of patterns 
is nicely revealed by the experimental work with pressure and axially 
loaded thin metal cylinders; Schnell (60) , Figure 5. 

Subject to extreme membrane penetration from high confinement 
there is in the early stages of loading a thin shell behavoir. The 
shell fails leaving surface irregularities in the form of lumps, 
bumps and pimples. It is important to note that preliminary 
calculations show the effect to be strictly superficial to the total 
loading-deformation behavoir, but as the loading progresses the surface 
irregularities grow and become exaggerated. In fact the growth is 
preferential with nonuniform geometry changes in that the "bulge" 
region irregularities show greatly accelerated growth. The local 
scale seems to predict and interact with the larger scales of 
deformation that are referred to as "bulge". 

These phenomena are shown on the length scales caricature; 

Figure 6, and are accompanied by a typical set of values as Figure 7, 
all related to the bead diameter. The first scale is a measure of 
the extent to which beads penetrate the confining membrane, with low 
penetration requireing a value 0(1) and is a measure of the angle 
through which the membrane turns in traversing a bead. Somewhat 
arbitrarily a membrane with t/d 0(0.1) will be termed thin and one 
with t/d 0(1) will be termed thick. 

As with a very uniform shell of homogeneous material, imperfections 
will greatly effect the load that marks catastrophe, (see Figure 10) 

Of course in the preparation of samples of material for a triaxial test 
leaves many fabric anomalies at the surface. These sites will 


XVII-12 


undoubtedly provide imperfections for buckling , and encourage premature 
failure. After failure these sites are marked with continual growth 
of vestigal Ivimps , ridges and depressions. 

The presumption is Lliat undei: cyclic loading the shell phencmenon 
would not reappear as the residual wrinkles would still be present. 
However they would without a doubt continue to grow monotonically 
under all cyclic loading being true catastrophe products and irrever- 
sible to the loading that produced them. Ironically, the postulate 
is that shellness is related to thin membranes; those with less internal 
bending stresses. The simple concept that hoop stress varies as 
the reciprocal of thickness must be viewed with askewance at high 
membrane strain. The extreme rise in modulus of elasticity with 
strain could cause thin membrane strains to be less than thick membrane 
strains; see Figure 11. 


Under the assumed sequence of membrane penetration proceeding 
axial load and the sxibsequent formation of an unbonded, friction 
linked shell; it is perhaps plausible that the shell would resemble 
the situation in a reinforced concrete shell. Seide(81) gives a 
broad survey, and using his equations as an estimation tool with the 
property values for latex rubber gives the results in Figure 9. The 
estimates are pres\imably lower bounds ; but as this shell does not 
have a distinct thickness or homogenuity of modulii there is no use 
to pursue this analog any further. The conclusion however is clear. 
In Figure 9 the "P" is the critical value of external pressure for 
buckling under that loading alone and the " fig " for axial stress 
alone; since both are typically present it might be of value to use 
Seide's combined loading figures, but to this level of approximation 
it is not helpful beyond the observation that critical axial loads 
would be reduced considerably. As an example consider the case of 
t+d = 3. 3mm.. the critical pressure of 0.01 psi must be related to the 
confinement minus the pore pressure ( or 1/3 tr ^ ) and the total 
axial load would be2i^Rt<J^ - 1.8 Lbf. The figure is intended to be 

schematic and does not show penetration which is assumed, t/d =0.1. 
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FIGURE 12 BUCKLED THIN SHELL 
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NUMERICAL APPROACH 


The ultimate source for the sequence of programs given the 
generic title MICROFEM is Chapter 24 ( by R.L. Taylor ) of the classic 
text by Professor Zienkiewicz, The Finite Element Method ^ McGraw (1982i 
The version produced in this work is a direct ancestor of several 
versions from the students of Professor Ko at the University of 
Colorado. Most recently DENISE. FTN by Garrett Denise, ( MSFC name). 

All are application oriented and have substantial differences, but 
the general solving subroutines are essentially unchanged. The most 
important feature being the Lade-Nelson constitutive model, added to 
the framework of the ultimate source. The next page is dedicated to 
listing, locating and typing the subroutines in this code. 


The program (MAIN ) drives 26 subroutines and 1 function. After 
an initial committment to nodalization, three subroutines are dedicated 
to checking size compatibility with the local environment: the current 
levels set are well below the capabilities of txie local PE 3250 on 
which the program is installed. The setup is done by a menu driven 
query in INPUT; this drives on command turning on of flags(l,7,8,9,10) 
direct entering of acceleration of gravity (11) , calling other 
Subroutxiies ( 2 , 3 , 4 , 5 , 6 , 12 ) and exit (13) ■ 


FORCE and AX I FORCE are used 


to commit the external loading; confining and axial respectively (6,12) 
Stress and strain at "Gauss points" is output from RESP01 which 
is called from several places in other subroutines via RESPON; the 
usual sequence is: LOAD, SOLVER, RESPON (RESP01 ) . 


The local name of MICROFEM is MGMFEM.FTN. It has been given 
a command file MGMFEM.CSS to drive it by its own name. Further a 
manual has been produced to provide some detail of the source and 
evolution for the MICROFEM family, details of operation, and coding 
flow to assist in modification. The manual includes; explicit startup 
instructions, hardcopy of the screens that appear during a run 
(ie: menus, prompts for input, mirrored array generation, and output), a 
subroutine listing with a brief description of utility, a subroutine 
level flow chart with ties to the data I/O points. There is also 
some command level instruction for generating a LOG file of an inter- 
active session, and the code for the autorun command file. 


The earlier versions of the code were run batch-style and setup 
from a data file. They have been retained locally as AXI.FTN, 

BI.FTN and DENISE. FTN. Also the BASIC program Capmodel has been 
converted from HP-BASIC to BASICA for the AT driving an HP plotter. 
There are 5 model stables in the MICROFEM code and 4 remain unused. 

One of the primary goals of the numerical phase of MGM analysis should 
be to populate those stables so that several types of constitution 
could be run with a single initialization. 
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LOCATION 


140 

• HI 
H2 

- H3 
H4 

• 343 

- 344 

■ 345 

■ 346 
485 

• 509 
516 

• 523 

• 524 

• 525 

• 612 

• . 679 

735 

801 

- 859 

■ 860 

■ 979 
980 

1010 

- 1011 

• 1128 

- 1129 

- 1170 

- 1171 

■ 1172 

- 1203 

- 1204 

- 1205 

- 1316 

- 1343 

- 1368 

- 1384 

- 1395 

- 1438 

- 1439 

- 162‘ 

- 1647 

- 1648 

- 1649 

- 1699 

- 1700 

- 1701 

- 1922 

- 1639 


SUBROOTINE LISTING: 


MGMFEM.FTN 


original page is 
OF POOR QUAUTf 


TYPE 


NAME 


ARGUEMENTS 


SOBRGUTINF: INITHI (TF tv 7n n y VI i:i 

1 NNP.HGL.NPlRT.NDin . NEN . NEH1 . NDOF , STRESS , Pfi . XKUR , XH , XNI . 
f .S1,S2.T1.T2.U1,U2.D1.C2,SK.BK.UBK.R,JDIA6. 

3 hSEC , NEC . NflC . IFfilL . TOL . DENSIT . NSTEP . NITER , 

4 PRFF.PRFC) 

SUBROUTINE rNCLOflnfTF . TV TP. n V F vi nt nnp nfi 
I mm ,MUir.NEN,NEN1 . NDOF . STRESS , Pfi . XKUR . XN , XNI . CL , PC , 

0 2?-^^,5MM2,T1,T2.U1.U2,01,C2,SK.BK.UBK,R,JDIHG,NSEQ, 

3 NEC . NfiD . IFfilL . TOL . NUniNC . DISP , PRFP . PRFC > 

SUBROUTINE CHKSTR a n 
SUBROUTINE IZERO (T.NUW) 

SjbKLIUTINE Rimj ( ft . NUN ) 

^■SitgRQUTINE INPUT(IE.3.IX.I0,X.F.IXL,XL.IDL.FL.NN.NNP.NEL. 

1 mft T . NO IN . NEN , NEN 1 . NDOF , Pfi . XKUR , XN . XNI , CL , PC . XN . ETft 1 . S 1 . 

.W2,C1 ,C2. TOL. NUNINC. DENSIT. NSTEP, NITER) 
lUBRCkTlNE. COORD < X . XL . NN . NNP . NOIN . IPRT , lERR ) 
subroutine FI TOPC IX . IXL . NEL . NEN . NENI . IPRT . IERR ) 

bUUKUUilNE FIXITY ( ID . IDL . NNP . NDOF . IPRT ) 

Subroutine . fl . id , nnp . ndof . iprt . ierr , nuninc > 

.X-XI NFI NFN NFN1 NflTM NnnP NNP 

N . IE . D . N6K . NGS . HNftT . IPRT , IERR . Pfi . XKUR , 
V,.^nA^!^i-^''-''^^’^'^-^'^^^*5’-S2.T1.T2.W1 .W2.Q1 ,02) 

SUBROUTINE NfiTOUN .IE.D.NGK.NGS NMftT TPPT Pfl yifliP VN vnt 
1 CL,PC,XN.ETfi1.S1.S2.T1,T2.U1.W2,Q1,02) 

.SUBROUTIN E PR0FIL( IX , ID . JDIfiG , NNP , NEL , NEN , NENI , NDOF . NST , 

1 NPu , NEu ) 

subroutine stiff ( IE . D . IX . X . ID , XL . SK . JDI fiG , BK . UBK , nnp , NEL , NN ftT . 

1 NEN NEN 1 . NO IN , NDOF . NED . NSED , NfiD , DSTRES , Pfi , XKUR , XN , XNI . CL , PC , 

2 .S2.T1 ,T2,W1 ,U2,Q1 .02, IFfilL . ITfiG, PRFP, PRFC) 

■ . B . IX . X . ID . XL . SK ■ JDIfiG ■ BK ■ UBK ■ NNP ■ NEL . 

1 HHfil .rthN.NLNI .NDIN , NDOF , NEC . NSED , NfiD , STRESS , Pfi , XKUR , XN , 

2 XNI ,CL.PC,XN,ETfi1 ,S1 .S2.T1 .T2,U1 ,U2,01 ,02, IFfilL, ITfiG, 
SUBROUTINE GfiUSOK I .1 TNT Sn Tfi LIG^ 

S UBROUTI N E SHfiP01 (S.T.XL.NDIN.NFN.ShP.R,nFT.n 

SK , LN , BK , UBK .JDIfiG , NSEO , NfiD , NED ) 
sub routine LOfiD( ID.F.R.NNP.NDOF.NFn) 

-SQLUEPI BK . UBK , R , JD I fiG . NED . NfiD . IFLfiG ) 

Subroutine DLfiDE< STR i , str2 . str3 . str 4 . Pfi , kur , n , ni , cl . pc . 

1 M.L |P1 ,S1 .SZ.TI ,T2,U1 .U2.01 .02. IFfilL. DEP,fiP.fiC,FP,FC) 
^I ' ePOUTINE gISECT P.P.O.FP Pft C P EP X TFI QR^ 
subroutine RESPONT IE. D. IX. X. ID. F. XL .U1 .R .NNP, NFI .NtlflT.NFN, 

1 NeNI .NDlli.NDOF.NSEQ . NEO . DDSTRE . STRESS . Pfi , XKUR , XN , XNI , 

2 CL.PC.XN.ETfil ,S1 ,S2.T1 ,T2.U1 .W2.01 .02, PRFP. PRFC, IFfilL) 
SUBROUTINE RESPOK IE.D.IX.X.F.XL .UL NNP.NFI NWflT.NFN.NFNI . 

1 NU i n . NUUh . UUS TRE . STRESS , Pfi , XKUR , XN . XNI , CL . PC , XN , ETfi 1 , 

'' S1.S2.T1.T2.W1 .U2.01 .02, PRFP, PRFC, IFfilL) 


0 

0 

s 

$ 


0 
0 
. « 
0 
0 
s 

0 

0 


$ 

0 


SUBRQUf INF PSTRFSTSffi PftTGV 

fGncTiiJn D0T<A.g7l 


* SETUP 

# GENERAL SOLVING 

$ MODEL RELATED 
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CONCLUSIONS & RECOMMENDATIONS 


The Jaumann Quandary has been laid out in some detail and put 
in perspective as a "state of the art" question in rational mechanics. 

The resolution of stress flux "choice" is pivotal to any confidence 
in numerical modelinq beyond the well established realms. A clear 
definition has been given for stress-power in the incremental sense. 
Apparently the impact of the stress-flux on this integral has not been 
realized before. The reoccurrence of the J.Q. has been shown at 
several places in the Scott (84) constitutive relation. 

The analytic path from optical displacement to the "true" strain (s) 
has been covered in some detail, but optical displacement is not yet 
resolved to numerical (FE) strain which will vary somewhat with specific 
algorithm. 

A close inspection should be made of the optimization process for 
t/d; the strategy to minimize menibrane thickness alone may be faulted . 
First: it is very important to control vestigal foirmations if surface 
strains are to be measured. 

Second: there is ultimately a thinness limit to lowered hoop stress 
established by the onset of increased modulus of elasticity. 
Third: the optical detection and selection process is much more robust 

thcin originally envisioned; and able to handle considerably 
higher t/d than 0.1. 

In light of the thin shell physical model presented above, the conclusion 
is that a thin membrane (t/d 0(.D) may encourage the production of 
vestigal lumps on the sample surface. These formations , while not 
changing the gross behavoir of the testing may cause considerable inter- 
ference with attempts to measure optical surface displacements. 

Hence the maximum t/d is encouraged which will still permit optical 
processes to remain error-free in following red-granule trajectories. 
Whether t/d can be made to approach 0(1) within the constraints of low 
confining pressure, minimum optical resolution and optical transparency 
remains to be seen. 

An improved version of "MICROFEM" has been produced; MGMFEM.FTN, 
which resides on the Fluid Dynamics Branch Perkin-Elmer 3250. It is 
fully interactive and can be flagged to mirror entry, output to screen 
or printer, reinput data, and so on. The data entry necessary for a 
run has been reduced by eliminating features not pertinent to the 
problem at hand; ie: RPM, other geometries, multiple material types, etc) , 
but the acceleration of gravity has been explicitly added as input. 

The following features still need to be modified; the code is still time 
insensitive, the loading is limited to stress initiation, the model 
stable has only Lade-Nelson, and many options for post-processing need 
to be implemented for a satisfactory analysis tool . 
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